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1. Introduction 



Recent development of the lattice QCD simulation enables us to access hadron properties 
with great accuracy. In fact, as a precise first principle calculation of QCD, the result of lattice 
QCD is now being used for examination of the fundamental theory, such as a unitarity triangle 
test of the standard model[[Ij]. Yet, there are various quantities remained unsettled so far. Among 
them, we focus on the isospin breaking effect on baryons, such as proton-neutron (p — n) mass 
difference. Actually, this mass difference is one of the most fundamental quantities in nuclear 
physics. For instance, this quantity governs the j3-decay of the neutron, i.e., the life time of the 
neutron. Note also that j8 -decay /electron capture are the basic ingredients in the understanding 
of the nucleosynthesis and the history of the universe: if the mass ordering between proton and 
neutron was opposite, the universe may not exist as it is now. In the laboratory-experiment, the 
isospin breaking effect on the baryon spectrum is observed not only for p — n mass splitting but 
also for splitting between other isomultiplets in octet/decuplet members. The charge symmetry 
breaking in the N — N interaction is also observed, which is caused by the isospin breaking of 
nucleons. 

The isospin breaking of hadrons originates in two ingredients: one is the up and down quark 
mass difference in QCD, and the other is the electromagnetic (EM) effect based on QED. In this 
sense, the determination of u, d quark mass through the study of isospin breaking corresponds 
to fixing the Yukawa coupling constant between u, d quark and the Higgs particle, which are the 
fundamental (and unknown a priori) parameters in the standard model. Moreover, u, d quark 
masses are particularly interesting from the viewpoint of the so-called strong CP problem. In fact, 
it is proposed that the existence of massless quark(s) can absorb the phase 6 of the vacuum, and 
thus resolve the strong CP problem. By studying the isospin breaking in hadron spectrum, we can 
examine whether or not such a scenario happens in our realistic world. On the other hand, the 
inclusion of QED in the QCD calculation is becoming the urgent task to develop the new physics 
search beyond the standard model. For instance, in the theoretical calculation of muon anomalous 
magnetic moment, it is know that hadronic light-by-light scattering contribution makes the large 
uncertainty, and it is proposed to resolve this problem using the QCD + QED simulation^. From 
this viewpoint, the study of isospin breaking on the spectrum is a suitable topic to build a firm 
foundation for the QCD + QED simulation. 

In the literature of lattice QCD, the first work has been done by Duncan et aZ.[|§], where 
they employed the quenched non-compact QED combined with quenched QCD with the Wilson 
fermion. In this pioneering work, however, both of the QCD quenching artifact and the artificial 
chiral symmetry breaking contamination are unavoidable. In this work, we eliminate these uncer- 
tainties by employing Nf = 2 dynamical domain-wall QCD configuration. Our study for the meson 
sector (such as 7t + — 7t°, K + — K° mass splitting) has been already reported in Ref.[|J], and we 
report the study of the EM effect on octet baryons in this proceeding. The effect of u, d quark mass 
difference will be reported elsewhere [Q]. There is another work|Q] on splitting between meson 
isomultiplets with improved action but still at the quenched level. Recently, the splitting of p — n is 
also studied|7]], while only the effect of u, d quark mass difference is considered there. 
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2. Formalism 

We study the EM effect on baryons through the mass difference between isomultiplets, such 
as m p —m n , m £ + — m^o, m^- — m^a, — w^-, m£+ + m^- — 2/m^o, m 5 - — m z o using the QCD 
+ QED simulation. In the QCD sector, we employ the Nf = 2 unquenched QCD gauge con- 
figuration D (;t) generated by the RBC collaboration^]. The domain- wall fermion action 
and the DBW2 gauge action is used with the parameters of V = L 3 x T = 16 3 x 32, L s = 12, 
M$ = 1.8, j8 = 0.8. The ensembles are generated for three different (u,d-) sea quark masses of 
fWsea = 0.02,0.03,0.04. The lattice unit is determined to be 1.691(53) GeV so as to reproduce p 
meson mass m p = 770MeV. The sea quark mass and the physical volume roughly correspond to 
m sea ~ \m s — m s (m s : strange quark mass) and L 3 ~ (1.9fm) 3 , respectively. We pick up about 200 
configurations from ~ 5000 trajectories available at each sea quark mass ensemble. The analysis 
is performed with bin size of 2, (i.e., bin size of 50 trajectory separation), in order to suppress the 
possible autocorrelation. 

In the QED sector, we employ a non-compact formulation of 
J^qed = jL.ji.v M^Ay 05 ^) — d v A^ ED (x)^j at the quenched level^]. In the generation, we first 

generate A^ ED in the momentum space under the Coulomb plus residual gauge fixing condition to- 
gether with a boundary condition for the constant modes. The configuration in the coordinate space 
is obtained by Fourier transformation. The advantage of this formulation is that the generation of 
A^ ED leads to just a Gaussian random number generation, and thus there is no autocorrelation be- 
tween the configuration even for arbitrary small coupling. In addition to that, we do not have to 
worry about the renormalization of the QED coupling constant, because the quenched QED in non- 
compact formulation is a free theory. Given the QED configuration described above, we obtain a 
U(l) link variable by U^ ED (x) = exp[— iA^ ED (x)], and then construct the QCD + QED configu- 
ration as U^ CD (x) x (£/^ ED (x))2, where Q = +2/3e, — l/3e for u and d quark, respectively. In 
order to study the QED charge dependence of the mass splitting, we use not only the physical QED 
charge, a em = e 2 /4n = o^), but also the charges of a em = (0), (except 
for a em = S^|L at m sea = 0.02). 

The mass of baryon B is measured using the two-point correlation function 
^BB{t) = Hx(Jb(x)Jb(Q)} with the use of positive parity projection. We use the operator Jb which 
has non-relativistic limit, for instance, = e a b c {uJ l Cj5db)u c as the proton operator. The operators 
for other octets and singlet are obtained by the SU(3) rotation. Although these procedures are valid 
for /7,7i,Z + ,r~,S~,S° baryons, additional treatment is necessary in the Q = 0,S = — 1 channel, 
i.e., for Z°,A8,Ai baryons. In fact, the mixing between Ag and Ai occurs because of the SU(3) 
breaking, and the mixing between £ ,Ag,Ai occurs because of the SU(2) breaking. Note here 
that E°( 1193) is massive than A(1116), experimentally. Therefore, in order to determine the mass 
splitting such as +m£- —Im^a, we have to extract the £° state as a first excited state in Q = 
0,S = —1 channel. For this purpose, we employ the so-called variational mefhod[|9|]. In practice, 
we calculate not only the diagonal correlation functions as = Ya(^Bf(x)JBi(0)) ,Bi = Bf £ 

{E°,A8, Ai} but also the off-diagonal correlation functions as HB f Bj(t),Bj ^ Bf G {lP, A&,A\}. By 
performing the diagonalization of the 3x3 correlation function matrix of n(? = ?o) 1 ■ n(0> we 
can extract not only the ground state A but also the first excited state £°. 
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3. Numerical results and discussions 

Using Nf = 2 unquenched domain-wall QCD configuration combined with the quenched QED 
configuration, we calculate the baryon two-point correlation functions for all octet and singlet 
baryons including the off-diagonal correlation functions in Q = 0,S = — 1 channel. By solving the 
inverse of the Dirac operator under the QCD + QED configuration, the EM effect is automatically 
taken into account for the valence quark. In order to suppress the contamination from excited states, 
we employ the wall source and point sink conelation function under the Coulomb gauge fixing. In 
this study, we consider the isospin breaking effect up to the first order. Therefore, in the evaluation 
of the EM effect, we can use the quark masses without the EM effect. Namely, we perform the 
lattice simulation only at the unitarity point, m U:Va \ ence = mj iVa i ence = m sea = 0.02,0.03,0.04 for u, 
d quarks, and we use m s = 0.0446 for s quark, which is determined from ^-input[|8|]. 

In the Q = 0,S = —1 channel, we use the variational method to extract the first excited state 
£ as well as the ground state A. As described previously, we first calculate 3x3 correlation 
matrix H(t) using the flavor bases of £°,A8,Ai, and then diagonalize Yl(t = to) -n(?), where 
the extracted eigenvalues correspond to the exponentially decaying correlation function for each 
ground/excited state. In the following analysis, we fix the arbitrary parameter to as to = 1. The 
dependence on to is discussed later. In the study of the mass splitting between £ triplets, we 
perform the similar procedure for £ + , £~ as well, i.e., we use Tlzz(t) /Utz(t = to) instead of simple 
njx(f) where £ = £ + ,£~. By employing this procedure, we can take the full advantage of the 
statistical correlation among £ triplets and obtain the reasonable signal. 

In order to study the QED charge dependence of the mass splitting, we evaluate the QED 
charges of a em = (0) , «|£ y , , , ML (except for a em = at m sea = 0.02). For each 

0£e m , we calculate the correlation function for not only e = +y/4na em but also e = —\/4na em and 
take the average between them. This corresponds to the use of the QED configuration of {A^ ED } 
— > {A^ ED , — A|J ED } with binning between A^ ET> and — A^ ED . In this procedure, &{e) contamination 
in the correlation function can be eliminated a priori. In fact, because the EM effect on physical 
observables appears only from (J(e 2 ), such &{e) contamination is nothing but a statistical noise. 
Practically, we find that this procedure actually improve the S/N drastically, and is essential for the 
study of the baryon EM splitting. 

In order to demonstrate the signal of the mass splitting, we analyze the ratio between the 
correlation functions of the isomultiplets. Considering the p — n mass difference for example, we 
can express the correlation function for proton and neutron as Yl pp (t) = A p exp[— m p t], Yl nn {t) = 
A„exp[— m„t], respectively, where X p is the proton (neutron) overlap constant between the 
state and the operator. Noting that (X p — X n ) and (m p — m n ) are Gie 1 ), we can write the ratio of the 
proton and neutron correlation functions as 

R p/n (t) = U pp (t)/U m (t) = 1 +2(X p -X n )/(X p + X n ) - {m p -m n )-t + 0{e A ). (3.1) 

Therefore, the slope of R p / n (t) in terms of the Euclidian time t directly corresponds to the p — n 
mass splitting. In Fig. [T], we plot R p j n (t) in terms of t. For each QED charge, we find a clear 
negative linear slope, which indicates m p > m n from the EM effect. 

In the practical calculation of the mass splitting, we perform the exponential fit for each baryon 
correlation function and evaluate the mass difference, where the statistical error is estimated by the 
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Figure 1: The ratio of the correlation function 
of proton and neutron R p / n {t) at m sea = 0.03 is 
plotted against t for each QED charge of e = 
(0),e p h y , 0.6, 0.85, 1.0. Negative linear slope corre- 
sponds to m p > m n from the EM effect. 
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Figure 2: The proton-neutron mass difference m p — 
m n from the EM effect at m sea = 0.03 is plotted 
against a em = <? 2 / (4ft). The result at aim is shown 
as the second point from the left. The solid line cor- 
responds to the best linear fit in terms of a em . 



jackknife method. In Fig. ||, we plot the p — n splitting determined at each QED charge a em . 
We observe that the splitting behaves linearly in terms of a em , and there is no indication of the 
appearance of higher order EM effect, ff(oc^ m ). Because we consider the EM effect up to &(a em ) 
as a framework, this observation guarantees that our procedure is self-consistent. 

Finally, we fit linearly the splitting in terms of a em , and determine the splitting at a em = aem y - 
We perform this procedure at each lattice simulation with quark mass of m = m va i enC e = ni&ea = 
0.02,0.03,0.04. Fig. [| shows the result for p — n mass splitting from the EM effect. We ob- 
serve the non-trivial EM effect on p — n for each m. Note that these are the first results obtained 
non-perturbatively using dynamical lattice simulation. It is interesting to see that the result at 
each m is roughly consistent with the model estimation using Cottingham formula, m p — m n = 
0.76(30)MeV[|§]. The splitting in the real world can be obtained by the chiral extrapolation of the 
lattice data. We, however, find that the result at m = 0.02 is afflicted with larger statistical fluctua- 
tion, and the reliable chiral extrapolation becomes difficult. In fact, the linear chiral extrapolation 
in terms of m leads to only zero-consistent result for p — n mass difference. In order to extract 
better signal, the statistical improvement is in progress. 

In Fig. I], we also show the result of the EM splitting for + m^- — 2m z o at each m = 
0.02, 0.03, 0.04. The splitting m z + + m z - — 2m z o is an interesting quantity from the viewpoint of 
discrimination of the two ingredients in isospin breaking, i.e., the EM effect and the u,d- quark mass 
difference effect. In fact, if we consider the SU(2) rotation of exp[— i7ia 2 /2] (i.e., u,d exchange), the 
£ triplets rotate as Z + — > L~ — > L + , E° — > E°. Therefore, in the splitting of m z + + m z - — 2m z a, 
there is no &{m u — mj) term and there exists only &(a em ) EM effect up to the first order of the 
isospin breaking effect. Under this consideration, the result from the EM effect can be directly 
compared with the experimental value of 1.5MeV, in principle. Unfortunately, because of the 
statistical fluctuation, we obtain only the zero-consistent result for m z + + m z - — 2m L o after the 
linear chiral extrapolation. The upper bound of the lattice result for this quantity is found to be 
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Figure 3: The proton-neutron mass difference from Figure 4: The mass difference + m^- — 2m z o 
the EM effect at physical QED charge = 1/137 from the EM effect at physical QED charge a|m Y 
plotted at each quark mass of m = m va i ence = m sea . plotted at each quark mass of m = »? V aience = '"sea- 



smaller than the experiental value, while this discrepancy may be attributed by the finite volume 
artifact[||]. 

Before closing this section, we comment on the systematic error in these results. In order to 
check the stability of the analysis, we examine the several alternative methods to extract the mass 
difference. For example, we perform the linear fit for the ratio of proton to neutron correlation 
functions (i.e., for the ratio shown in FigJI]) in terms of t, instead of exponential fit of each correla- 
tion function. We find that the result is consistent with each other and confirm that the results are 
reliable. In the variational method, we choose several ?o instead of to = 1 choice when diagonaliz- 
ing Yl(t = to)~ l • 11(7), and check the dependence on to- It is found that the results are insensitive 
to to, and we confirm that our variational procedure is stable. Yet, the current lattice results are 
afflicted by the statistical noise, particularly in light quark mass sector, and it is difficult to perform 
the definite chiral extrapolation. Further calculation is desirable to achieve the better statistics, 
which is actually our ongoing work. Although the uncertainty in the result extracted from the 
current lattice setup has been evaluated as describe above, the most troublesome artifact remained 
is the finite volume artifact, because the QED interaction is a long-range interaction. In fact, the 
model-based calculation^] suggests that such artifact is sometimes comparable to the results of 
the lattice simulation. At this moment, we cannot evaluate the finite volume artifact without using 
model calculations, and the results given above would receive some modifications. The explicit 
calculation of the finite volume artifact is planned with the use of Nf ■ = 2 + 1 configurations gener- 
ated by the RBC-UKQCD collaboration [JTTJ], with which the configurations with different volumes 
are available. 

4. Summary and outlook 

We have investigated the electromagnetic (EM) effect on octet baryon spectroscopy. By em- 
ploying Nf = 2 dynamical domain-wall QCD configuration combined with non-compact quenched 
QED configuration, the u, d sea quark effect in QCD has been incorporated. The mass split- 
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tings between isomultiplets (p,n), (E ,E~) have been studied by evaluating the two- 

point correlation function, where the variational method has been adopted in the Q = 0, S = — 1 
channel in order to extract the E° state as a first excited state in this channel. In order to study 
the QED charge dependence of the mass splitting, we have chosen the QED charges of a em = 
(0), 0(£m , 4^ > ^jjr i ^4 J ■ We have found that it is essential to calculate both of e = ±\/47ia em 
for each a em , in order to cancel the G{e] contamination and to achieve the reasonable S/N in 
the baryon EM mass splitting. By fitting the mass splitting linearly in terms of a e m> we have 
obtained the first result from lattice dynamical simulation for the baryon EM mass splitting at 
each m = 0.02,0.03,0.04. The investigation of the effect of u, d quark mass difference is also 
in progress^]. There remains uncertainty originates from the finite volume artifact. In order to 
investigate it explicitly, we are planning to perform the analysis using the Nf = 2 + 1 dynamical 
domain- wall configuration generated by the RBC-UKQCD collaboration[[II]]. There, we can also 
eliminate the quenching artifact for the strange quark as well. In future, the inclusion of the dy- 



namical QED effect is interesting to investigate[]12|], in which we expect the definite calculation of 



the isospin breaking effect is possible with all the uncertainties in the simulation under control. 
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